Functional renormalization for spontaneous symmetry breaking in the Hubbard model 



S. FriederichEl H. C. Krahl, and C. Wetterich 

Institut fur Theoretische Physik, Universitdt Heidelberg, Philosophenweg 16, D-69 120 Heidelberg, Germany 

The phases with spontaneously broken symmetries corresponding to antiferromagnetic and J-wave super- 
conducting order in the two-dimensional t - -Hubbard model are investigated by means of the functional 
renormalization group. The introduction of composite boson fields in the magnetic, charge density and super- 
conducting channels allows an efficient parametrization of the four-fermion vertex and the study of regimes 
where either the antiferromagnetic or superconducting order parameter, or both, are nonzero. We compute the 
phase diagram and the temperature dependence of the order parameter below the critical temperature, where 
antiferromagnetic and superconductiving order show a tendency of mutual exclusion. 
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I. INTRODUCTION 

The two-dimensional Hubbard model on a square lat- 
tice is widely believed to hold a key role for the understand- 
ing of high-temperature superconductivity in the cuprates. In 
analogy to the cuprate phase diagram it seems to exhibit an- 
tiferromagnetic and J-wave superconducting order in close 
vicinity |4-20]; for a systematic overview, see Ref. [211. The 
tendency toward J- wave superconductivity was already pre- 
dicted by some strikingly simple scaling approaches ll22l - [24l] . 
On a higher level of technical sophistication, the fermionic 
functional renormalization group approach [25-32] has been 
of great help to analyze in detail the competition of different 
types of instabilities. 

The main result of the present paper concerns the phase 
diagram in the low temperature regime, where spontaneous 
symmetry breaking occurs either in the antiferromagnetic or 
in the J- wave superconducting channel. The symmetry bro- 
ken regimes are difficult to access for fermionic functional 
renormalization group studies such as Renormal- 
ized mean field investigations lHH [sj], based on a combi- 
nation of the fermionic functional renormalization group ap- 
proach with a mean field treatment, have already been able to 
study the mutual influence of the order parameters for antifer- 
romagnetism and J-wave superconductivity. The functional 
renormalization group approach presented in this paper makes 
it possible to investigate this problem by entering the spon- 
taneously broken regimes in a partially bosonized language 
where the diflTerent types of bosonic fields that are introduced 
correspond to the diflferent types of possible order of the sys- 
tem. 

We find a region of electron fillings with competition be- 
tween antiferromagnetic and superconducting order. The two 
phenomena show a strong tendency of mutual exclusion. Nev- 
ertheless, we find regions where the two types of local order 
coexist, even though this coexistence may not be maintained 
on the length scales of global order. We also find near the 
van Hove filling a considerable range of temperatures with 
local but not global antiferromagnetic order. For smaller elec- 
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FIG. 1: Phase diagram for U/t = 3 and f It = -0.1. The (blue) solid 
region for smaller values of denotes the antiferromagnetic phase, 
the (red) solid region for larger values of the J- wave superconduct- 
ing phase. The green dashed line and the black dotted line indicate 
the pseudocritical temperatures below which local order sets in for 
antiferromagnetism and J- wave superconductivity, respectively. The 
pseudocritical line for antiferromagnetism ends at -fi/t ^ 0.79. The 
region below T^nm = 4 • 10~^^ has been left blank since calculations 
were only done for higher temperatures. The range of the chemical 
potential fi shown here corresponds to electron filling between 0.9 
(for -jd/t ^ 0.4) and 0.7 (for -ji/t ^ 0.9) per site. (We have no 
precise estimate for (rie) as a functiona of fi for all T.) 



tron filling the critical temperatures for the onset of local and 
global superconducting order almost coincide. An overview 
of our findings in form of the yU-T-phase diagram (for next- to- 
nearest neighbor hopping f /t = -0.1 and Hubbard repulsion 
U/t = 3) is shown in Fig. [T] 

The approach presented in this work focuses on the low 
temperature behavior and brings together earlier attempts 
1.35-411 to perform a functional renormalization group anal- 
ysis of the two-dimensional Hubbard model based on par- 
tial bosonization (or Hubbard- Stratonovich transformation) 
JH, liH. In analogy to the parametrization methods for 
the fermionic four-point vertex proposed and developed in 
iH. lisll this method makes it possible to treat the complex 
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momentum dependence of this function in an efficient, simpli- 
fied way, involving only a comparatively small number of cou- 
pled flow equations. The fermionic four-point vertex, which 
is a scale-dependent function of three independent momenta, 
is decomposed in terms of bosonic propagators and Yukawa 
couplings, which are each functions of only one variable. 

The main advantage of the method used in the present 
work consists of the possibility of following the renormal- 
ization group flow into regimes where one or more symme- 
tries of the Hubbard Hamiltonian are spontaneously broken. 
(For other renormalization group studies of symmetry bro- 
ken phases in similar models see ll4^ - [48ll .) At a certain scale 
kssB of the renormalization flow, the momentum-dependent 
fermionic four-point vertex may diverge, and this signals the 
onset of local collective order. In order to extend the renormal- 
ization group treatment to the locally ordered regimes, com- 
posite degrees of freedom such as magnons and Cooper pairs 
are made explicit in terms of composite bosons. These corre- 
spond to difl'erent types of collective order the system might 
exhibit. A nonzero expectation value of the magnon field, for 
instance, signals the presence of some form of magnetic or- 
der, and a nonzero value of some Cooper pair field signals 
superconducting order. Diff'erent Cooper pair fields are dis- 
tinct due to diff'erent symmetries of the order parameters they 
correspond to. The language of partial bosonization, where 
the diff'erent types of bosons are taken into account explicitly, 
is therefore the right tool to investigate the regimes exhibiting 
diff'erent forms of collective order. A particular advantage of 
the present approach, which combines functional renormal- 
ization and partial bosonization, is that it allows to investigate 
the possible coexistence of diff'erent types of order in the same 
range of parameters. 

Earlier renormalization group studies using the framework 
of partial bosonization (see Refs. [s^ [stI [39l incorpo- 
rate only a comparatively small number of bosonic fields and 
therefore obtain a poorer resolution of the four-fermion ver- 
tex. In addition, they suff'er from the so-called mean-field am- 
biguity 1 49] which refers to the fact that by the use of the Fierz 
relations the microscopic Hubbard action can be mapped by 
diff'erent Hubbard- Stratonovich transformations into various 
equivalent descriptions involving diff'erent bosonic fields. In 
the presence of approximations the final results eventually de- 
pend on the specific choice of the decomposition. As argued 
in Ref. 41, this shortcoming can be avoided in the present 
approach by explicitly keeping the original Hubbard repul- 
sion U in the truncation. Contributions to the four-fermion 
vertex that emerge during the flow are attributed to the dif- 
ferent bosonic channels, depending on their momentum struc- 
ture and therefore in an essentially unbiased way. A remain- 
ing bias concerns the choice of bosons taken into account. We 
choose to include four difl'erent bosons corresponding to mag- 
netic, charge density, 5-- wave and J- wave superconducting or- 
der. This choice can be motivated by the structure of the one- 
particle-irreducible (IPI) diagrams from which the contribu- 
tions to the flow of the fermionic four-point vertex are derived 

EH]. 

Our truncation for the flowing action includes, for the sym- 
metric regime, the parametrization of the four-fermion ver- 



tex already used in Ref. 41. Furthermore it contains quartic 
bosonic couplings for the antiferromagnetic and J-wave su- 
perconducting bosons as well as self-energy corrections for 
the fermions in terms of a frequency-dependent wave func- 
tion renormalization. We observe symmetry breaking within 
the studied parameter region only in the antiferromagnetic and 
J- wave superconducting channels. In these symmetry broken 
regimes, the p- and ^-bosons are dropped from the truncation 
for renormalization scales k < kssB, and we focus on the 
fermionic and bosonic contributions to the order parameters 
and quartic couplings of the antiferromagnetic and J- wave su- 
perconducting bosons. 

II. METHOD AND APPROXIMATION 

The starting point of our treatment is the exact flow equa- 
tion for the eff'ective average action or flowing action [51]: 

dkTk = ^STr (rf + Rk)"' ^kRk = ^STr^k ln(rf + Rk) . 

(1) 

The dependence on the renormalization scale k is introduced 
by adding a regulator Rk to the full inverse propagator F^^^ . In 
Eq. ([T]) STr denotes a supertrace, which sums over momenta, 
frequencies, and internal indices, containing a minus sign for 
fermions, while dk = (dkRk)-^ is the scale derivative acting 
only on the scale dependence introduced by the regulator Rk. 
The Hamiltonian of the system under considerations is taken 
into account by the initial condition Fy^^A = 5" of the renormal- 
ization flow, where A denotes some very large UV scale and S 
is the microscopic action in a functional integral formulation 
of the Hubbard model. For the Hubbard model, this action is 
given by 

S = Yj^\Q)[iojQ^^Q]HQ) (2) 

Q 

xS(Ki -^2 + ^3-^4) , 

where 

HQ) = {h(Q)^HQ)f (3) 

are Grassmann fields describing electrons on a square lattice. 
The next-neighbor and next-to-nearest-neighbor hopping pa- 
rameters t and f are reflected in 

^(q) = -jd- 2^(cos qx + cos qy) - 4f cos q^ cos qy . (4) 

Here, as well as in all that follows, we employ a compact 
notation with Q = (ojn = 27inT, q) for bosonic and Q = (ojn = 
(2n -h l)7iT, q) for fermionic fields, and 

Q n--oo ^ ^ ^ 

^ -7: 

m - Q') = r-'5„,„.(2;r)V2>(q - q') . (5) 
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The components of the momentum q are measured in units of 
the inverse lattice distance a"^ The discreteness of the square 
lattice is reflected by the 27r-periodicity of the momenta q. 

In the limit ^ ^ all fluctuations are included and the 
flowing action Tk equals the full efl'ective action T = Tk^o, 
which is the generating functional of the one-particle irre- 
ducible (IPI) vertex functions. For ^ > the bosonic fluc- 
tuations with momenta |q| < ^ and the fermionic fluctuations 
with momenta |q - q^l < where qp is the Fermi-momentum 
vector which is closest to q, are not yet included. 

Although Eq. ([T]) is an exact flow equation, it can only 
be solved approximately. In particular, a truncation has to 
be specified for the flowing action, indicating which of the 
(infinitely many) IPI vertex functions are actually taken into 
account. We use diff'erent truncations for the disordered sym- 
metric regime (SYM) and the spontaneously broken regimes 
(SSB) where one of the bosonic fields has a nonzero expecta- 
tion value. In what follows we denote the regime where both 
the order parameter ao for antiferromagnetism and for d- 
wave superconductivity are nonzero by SSB ad. The regimes 
where only either ao or is nonzero, the other one being zero, 
are denoted by SSB a and SSBd, respectively. The nonzero 
expectation values ao or in the SSB regimes indicate local 
order. Since we are dealing with a two-dimensional model, the 
order parameters ao and must become zero in the thermo- 
dynamic limit for ^ i-^ 0, in accordance with the Mermin- 
Wagner theorem. They may, however, remain nonzero for 
k < /"\ with / the size of a typical experimental probe, sig- 
naling the appearance of long-range, that is "global", (for in- 
stance magnetic or superconducting) order on macroscopic 
length scales. 

Our ansatz for the flowing action includes terms for the 
electrons, for the bosons in the magnetic, charge, ^-wave and 
J- wave superconducting channels, and for interactions be- 
tween fermions and bosons: 

^klx] = ^F,k + ^Fm,k + ^Fp,k + ^Fs,k + ^Fd,k (6) 

+^a,k + ^p,k + ^s,k + ^d,k + ^ UB,k(^,p, d) . 

X 

The collective field ;^ = (a,p, s, s*, d, d*, i//, i//*) includes both 
fermion fields i//, i//* and boson fields a,p, s, s*, d, d*. The "an- 
tiferromagnetic" boson field a( Q) is related to the "magnetic" 
boson field m(2) used in Ref. |4ll by a shift in the momen- 
tum variable with respect to the antiferromagnetic wave vector 
n = (0,7r,7r), 

a(G) = m(e + n). (7) 

The purely fermionic part Tf (the dependence on the scale 
k is always implicit in what follows) of the flowing ac- 
tion consists of a fermion kinetic term F/rkin, a momentum- 
independent four-fermion term F^, and the momentum- 
dependent four-fermion terms F^, F^, F^, F^: 

Tf = F^kin +F^ +F^ + F^ +F^ +F^ . (8) 

The fermionic kinetic term is essentially left unchanged with 
respect to the original Hubbard Hamiltonian dJ]), apart from 



the fact that a fermionic wave function renormalization is in- 
cluded, which depends on the Matsubara frequency, 

rFy.n = J]i/^\Q)PF(Q)HQ), (9) 

Q 

where 



PF(Q) = ZF(0JQ){iC0Q^m) . (10) 

On initial scale ^ = A we set Zf(ojq) = 1 for all frequencies 
in order to equal the kinetic term of the microscopic Hub- 
bard action, see Eq. The flow of Zf(cl>q) is neglected for 
all frequencies except for the two lowest Matsubara modes 
ojQ = ±7iT. The computation of the scale-dependent quantity 
Zf(±7iT) is described in the following section. Self-energy 
corrections to the dependence of Pf(Q) on spatial momentum 
are omitted. According to llsol] the influence of self-energy 
corrections due to the frequency dependence of Pf(Q) seems 
to be more important. 

As motivated in the introduction, the momentum- 
independent part of the four-fermion coupling, which at ^ = A 
is identical to the Hubbard interaction U, remains unmodified 
during the flow. The corresponding part of the efl'ective action 
reads in our truncation 

^F = \ Yj U6{Ki-K2 + K^-K,) 

X [ilf\Ki)ilf{K2)] [ilf\K^)ilf{K,)] . (11) 



We focus on the momentum and spin dependence of the 
fermionic four-point function Af(Ki,K2, K3, K4), which, due 
to energy-momentum conservation, is a function of three in- 
dependent momenta (e.g., K4 = Ki + We decompose 
this vertex into a sum of four functions /l^(2), /1^(2), Ap{Q) 
and A^p(Q), each depending on only one particular combina- 
tion of the Ki. The chosen decomposition of the fermionic 
four-point function is inspired by the singular frequency and 
momentum structure of the leading contributions during the 
renormalization flow, see Eqs. (9)-(12) in Ref. 41 for precise 
definitions. In our approach, these functions are described by 
the exchange of the four diff'erent bosons a, p, s and d. 

Practically, this is achieved by the technique of flowing 
bosonization |[52l - [54ll . which was adapted to our purposes in 
Refs. [3^ llil. The basic idea is to introduce scale-dependent 
bosonic fields in order to transform the momentum-dependent 
four-fermion vertexinto Yukawa-type interactions between the 
fermions and bosons.In this way we keep the terms F^, F^, 
F^ and F^ at zero during the flow and replace their efl'ects by 
flowing Yukawa interactions between fermions and bosons. 
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These interaction terms read in our truncation 

Tfa = '^-(^ «(^) • l'/^^(Q)'^W)] 



K,Q,Q' 



6(K-Q + Q' +n), 
r^-p = - E '^P^^^ P^^^ W\Q)m')] S(K-Q+ Q') , 



K,Q,Q' 



(12) 



K,Q,Q' 



-s{K) [ilfHQ)er(Q')]) 5{K-Q- Q') , 
Tpd = hdiK)fd {{Q- Q')l2){d\K)[V{Q)eilf{Q')] 



K,Q,Q' 



-d(K) [i/f\Q)ei/,*(Q')]) 6{K-Q- Q') , 



where cr = (cr^ , cr^ , cr^Y is the vector of the PauH matrices 
and 6 = icr^. The Yukawa-coupHngs ha(Q),hp(Q),hs(Q), and 
hd(Q) are running coupHngs which vanish on the initial scale 
^ = A of the renormalization flow since there is no non-trivial 
momentum dependence in the initial four-fermion term in the 
Hubbard action Eq. Q. Note the presence of the d-wave form 
factor 

fd(Q) = /j(q) = ^ (cos(^J - cos(^3;)) (13) 

in the second-to-last line of Eq. ([T2l) , which is kept fixed on 
all scales. 

The purely bosonic part of our truncation for the eff'ec- 
tive average action consists of the bosonic kinetic terms to- 
gether with the bosonic eff'ective potential. The kinetic terms 
of the bosons are defined as the momentum-dependent pieces 
Pi(Q) of the inverse bosonic propagators. The inverse propa- 
gator of, for instance, the antiferromagnetic boson is given by 
Pa(Q) = Pa(Q)-^^l, where fhl is its minimal value and Pa(Q) 
the (strictly positive) kinetic term which we parametrize as 



Pa(Q)=ZaOjl^AaF(q). 

In this equation we employ for F(q) 

Dl • [q]2 



Dl + [q]2 ' 



(14) 



(15) 



if commensurate antiferromagnetic fluctuations dominate. 
Here [q]^ is defined as [q]^ = -\- for qxj e [-7r,7i] and 
continued periodically otherwise. If incommensurate antifer- 
romagnetic fluctuations dominate, we use 



Fi(q,q) 



Dl^F(q,q) 



(16) 



where the momentum dependence is quartic in momentum 
and explicitly includes the incommensurability q: 



m q) = ^((^' - [q]')' + ^VqxfVqyf) • 



(17) 



The shape coefficient Da used in Eqs. (fTSl) and ([T6l) is defined 



as 



Da= ^{Pa(0,7T,7T)-Pa(0,q,0)). 



(18) 



For the prescriptions used in the computation of and in 
the symmetric regime and the parametrizations of the kinetic 
terms of the other bosons see Eqs. (B6) - (B8) in Ref. 41. No 
linear frequency term is included in the kinetic terms of the 
superconducting bosons. 

The contributions to the eff'ective average action where the 
bosonic kinetic terms appear are 



= ^y.^^(-Q)Pa(QMQ), 

Q 

^ Q 

= 2^*^2)/^^^2)^^2), 

Q 

= Yj^\Q)Pd{Q)d{Q). 



(19) 

(20) 
(21) 
(22) 



One can reconstruct the momentum-dependent four-fermion 
interactions by solving the field equation for the bosons / 
as a functional of fermionic variables (as derived by variation 
of F with respect to the field for the boson /) and reinserting 
this functional into F. Our results for the fermionic four-point 
function are one-loop exact in the sense that the scale deriva- 
tives of all contributions up to second order in the Hubbard 
interaction U are taken into account by our truncation, includ- 
ing their full dependence on spatial momentum. 

We may summarize that in one-loop order the complicated 
spin and momentum dependence of the fermionic four-point 
function, as it emerges during the renormalization flow, is 
completely expressed by the bosonic propagators and Yukawa 
couplings connecting the fermions to the difl'erent bosons. We 
expect that also beyond one-loop order the dominant features 
of the momentum dependence of Ap are reasonably well re- 
produced by the solution of the flow equations in our trunca- 
tion. 

We also include in our truncation a local efl'ective potential 
U b(^,P, s, d) (not to be confused with the Hubbard interaction 
U). Here we make an expansion in powers of fields a, p, s, 
and d up to second order in p and s and up to the fourth order 
in a and d. This expansion has its limitations. For instance, 
it cannot describe first order transitions between two diff'er- 
ent phases with the same symmetries. For our purposes the 
polynomial expansion is expected to work reasonably well. 
Spontaneous symmetry breaking in the antiferromagnetic or 
superconducting channels can be described by a minimum of 
Ub away from the origin in a-J-space. In case of a second 
order phase transition this means that the terms quadratic in 
the fields a and J, J*, evaluated at a macroscopic scale kph, 
turn negative for temperatures below the critical temperature 
T < Tc and vanish for T = Tc. We denote the quartic cou- 
pling in the antiferromagnetic channel by Aa, the coupling in 
the J- wave superconducting channel by Ad, and the coupling 
between these two channels by Aad- In the symmetric regime 
S YM we expand the eff'ective potential around the zero value 
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of the fields: 

Y,UB{2i,p,s,d) = ^i(m2a^(-e)a(e) + m>(-e)p(e)) 

X Q ^ 

+m]s\Q)s{Q) + mld\Q)d{Q) 

+ \ Yj ^(21 + 22 + 23 + 24) 

x(^a(ei, 22^23, 24) (23) 
+^^^(2i, 22)^(23, 24) 
+2^^^(21,22)^(23,24)) , 

where we have defined the quantities o^(2i, 22) = \^{Qi) ' 
a(22) and d{Qu Qi) = d\Qi)d{Q2) (which has to be distin- 
guished from the Dirac delta- function by the number of argu- 
ments). 

In the spontaneously broken regime SSBad the minimum 
of the eff'ective potential occurs at nonzero values of the fields 
a and d. In this case, we neglect the p- and ^--bosons in our 
truncation and expand around the minimum of the eff'ective 
potential at (ofo, ^o)^ 

^ UB{2i,d) = i Yj ^(2i + 22 + 23 + 24) 

X Qi,Q2,Q3,Q4 

{AaHQuQ2)-aoS(Qi)S(Q2)} 

+M^(Qu Qi) - ^0^(20^(22)} (24) 

X{^(23,24)-M(23)^(24)} 

+2^^{a(2i, 22) - ocMi)mi)\ 

X{c5(23, 24) -^0^(23)^(24)}). 

In the regimes SSBa and SSBd, where only either or is 
nonzero, the mass term for the boson with zero order parame- 
ter is kept in the truncation for the eff'ective potential. 

The parametrization we use for the frequency- and 
momentum-dependence of the bosonic propagators and the 
Yukawa couplings can be found in Appendix B of Ref. lU 
The sole diff'erence between the truncation used here for the 
symmetric regime SYM and in Ref. 41 is that nonzero quartic 
bosonic couplings Xa, Ad and Aad and a fermionic wave func- 
tion renormalization factor ZpinT) are taken into account in 
the present work. The presence of the quartic bosonic cou- 
plings is crucial for the flow in the symmetry broken regimes. 
In a purely fermionic language they correspond to vertices 
with eight fermions. 



III. INITIAL CONDITIONS AND REGULATORS 

As the microscopic scale k = A goes to infinity, the flow- 
ing action must be equivalent to the microscopic action of 
the Hubbard model, so the initial value of the four-fermion 
coupling must correspond to the Hubbard interaction U. The 
bosonic fields decouple completely at this scale, where the 



initial values of the Yukawa couplings are given by 

ha\A = hp\A = hs\A = hd\A = . (25) 

In practice, we choose a finite but very large A, which is a 
very good approximation. 

For the bosonic mass terms we take fhf^ = f and P/,a = 0. 
The choice m^^ = f amounts to an arbitrary choice for the 
normalization of the bosonic fields, which are introduced as 
redundant auxiliary fields at the scale k = A, where they do 
not couple to the electrons. Of course, this changes during the 
flow, where the bosons are transformed into dynamical com- 
posite degrees of freedom, with nonzero Yukawa couplings 
and a nontrivial momentum dependence of their propagators. 
The quartic bosonic couplings vanish on initial scale k = A. 

In addition to the truncation of the efl'ective average ac- 
tion, regulator functions for both fermions and bosons have 
to be specified. We use "optimized cutoff's" (551 for both 
fermions and bosons. The regulator function for fermions is 
given by 

Rk (2) = sgn(^(q)) (k - l^(q)l) 0(k - |^(q)|) , (26) 
while the regulator functions for the real bosons are given by 

K^'iQ) = Aalp • {k^lt^ - Fe//(q, q))m^/t^ - Fc/i(q, q)) (27) 

allowing for an incommensurability q with Fcn as defined in 
Appendix B of Ref. 41 (with Aa = Ajn and an additional H- 
shift for the <2-boson). Regulator functions for the Cooper-pair 
bosons are of the same form, but no incommensurability needs 
to be accounted for in these cases. 



IV. FUNCTIONAL RENORMALIZATION FOR THE 
SYMMETRIC REGIME 

The ffow equations for the couplings follow from projec- 
tion of the exact ffow equation for the eff'ective average action 
onto the various diff'erent monomials of ffelds. The right hand 
sides of these ffow equations are given by the IPI diagrams 
having an appropriate number of external lines, including a 
scale derivative dk = (dkRk)-^^ acting only on the IR regula- 
tor Rk. Diagrams contributing to the ffow of boson propaga- 
tors and Yukawa couplings which do not include any quartic 
bosonic couplings have been discussed in Ref. 41. Here we 
focus our discussion on diagrams which contribute to the ffow 
of the quartic bosonic couplings Aa, Ad and Aad, and on how 
these couplings aff'ect the ffow of the bosonic mass terms and 
Yukawa couplings. We neglect the quartic couplings for the ^- 
and p-bosons since the corresponding channels do not exhibit 
critical behavior in the parameter regimes studied. 
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FIG. 2: Diagrams involving the quartic bosonic coupHngs Aa, Ad and 
Aad which contribute to the flow of the antiferromagnetic and J- wave 
superconducting propagators in SYM. Scale derivatives of the dia- 
grams in the first line contribute to the flow of the propagator for the 
antiferromagnetic spin waves, those in the second line to the one for 
J-wave superconductivity. Wiggly lines denote antiferromagnetic, 
dashed lines superconducting bosons. 
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FIG. 3: Flow of the antiferromagnetic mass term ml for U/t = 3, 
t' It = -0.1, iJLit = -0.77 and T/t = 0.0215. The inset shows a 
detail of the flow where reaches its minimal value, followed by 
an increase due to the bosonic contributions in the second line of Eq. 
(ESI. 



A. Bosonic mass terms 

The flow of the antiferromagnetic and J-wave supercon- 
ducting mass terms is given by the following two equations, 

kdkfhl = 2hl(0) y kdk- (28) 

-\ kd I \ 

and 

kdkml = -4hl{0) y kdk , ^'^'"^f (29) 

Zkd {2 + ^ 1 

The first and second lines of these equations correspond to 
the fermionic and bosonic loop contributions, respectively. In 
Fig. [2] we show a graphical representation of the contribution 
from bosons, while the fermionic diagrams can be found in 
Fig. 2 of Ref. 41. The momentum vector Q, which appears 
in the denominator of the fermionic contribution to kOkth^ in 
Eq. (|28l) , accounts for the dominance of incommensurate over 
commensurate antiferromagnetic fluctuations in a wide range 
of parameters. In this case the minimum of the inverse an- 
tiferromagnetic propagator Pa{Q) + no longer occurs for 
2 = 0. Rather, there exist four discrete minima at vectors 
±Qx = ±(0, 0) and ±Qy = ±(0, 0, q), either of which can be 
used as the vector Q in Eq. (l28l) . For a detailed description 
of how incommensurate antiferromagnetism is treated within 
the present approach see Ref. HO- 

In the symmetric regime, the fermionic contributions de- 
crease the bosonic mass terms during the flow, whereas the 
bosonic contributions, proportional to the quartic couplings 
Aa, Ad and Aad^ tend to increase them for positive Aa, Ad and 
Aad' The closer the mass terms approach zero, the more im- 
portant the bosonic fluctuations become. Once the a- or d- 



boson mass term becomes close to zero, bosonic fluctuations 
become relevant, and the terms in the second lines of Eqs. ([28l) 
and ([29l) may prevent it from actually reaching zero. Fig. [3] 
shows this by means of an example where the bosonic contri- 
bution proportional to Aa inverts the direction of the flow of 
the mass term so that it remains nonzero for ^ ^ 0. 

Whenever a bosonic mass term becomes zero during the 
flow, we change our description of the eflTective potential from 
the form of Eq. (l23l) to that of Eq. (l24b or the corresponding 
versions for SSBa and SSBd. A negative quadratic term in the 
eflTective potential indicates local order, since at a given coarse 
graining scale k the eflTective average action evaluated at con- 
stant field has a minimum for a nonzero value of the boson 
field. The largest temperature where at fixed values of U,t\ii 
one of the mass terms vanishes during the flow is called the 
pseudocritical temperature Tpc. It can also be described as the 
largest temperature where short-range order sets in. At this 
temperature the eflTective momentum-dependent four-fermion 
coupling diverges in the channel where hits zero, as seen 
in the purely fermionic flow studies like Ref. |28|. However, 
the local order does not necessarily lead to long-range order, 
since the tendency toward order may be countered by long 
range bosonic fluctuations. 

If the order persists for k reaching a macroscopic scale, the 
model exhibits eflTectively spontaneous symmetry breaking, 
associated in our model to (either commensurate or incom- 
mensurate) antiferromagnetism or J- wave superconductivity. 
The true critical temperature Tc is defined as the largest tem- 
perature for which local order persists up to some physical 
scale kph corresponding to the inverse size of a macroscopic 
sample (see Refs. |3^|38|. We choose here kph = {lcm)~^ ^ 
lO'^t. In order to determine the true critical temperature for 
either a- or J-type of order, it is therefore necessary to switch 
to the truncation in which either ao or (or both) are nonzero. 
Already a quick inspection of the phase diagram forU/t = 3 
and f /t = -0.1 (Fig. [B reveals the importance of the flow 
in the spontaneously broken regimes. The pseudocritical tem- 
perature Tpc diflTers substantially from the critical temperature 



-5 5 

-ln(k/t) 

1.5 I ^ ^ r 




-5 
-ln(k/t) 

10 I ^ ^ 




-5 5 

-ln(k/t) 



FIG. 4: Flow of the bosonic mass terms ml, m^, and (upper 
panel), the Yukawa couplings ha(0) , Zip (11), hs(0) and hd(0) (mid- 
dle panel). The lower panel shows a logarithmic plot of the ef- 
fective fermionic four-point couplings Ap^ where Af^a = ^1(0)/ ml, 
Af,p = hl(U)/ml, Af,s = hl(0)/ml and Af,d = hl(0)/ml. The lines for 
all three panels are (green, dashed) for the antiferromagnetic boson, 
(red, solid) for the J- wave superconducting boson, (blue, dotted) for 
the charge density wave boson, and (magenta, dashed-dotted) for the 
5-wave superconducting boson. All dimensionful quantities are in 
units of t. Parameters chosen are U/t = 3, f /t = -0.1, fi/t = -0.6 
and T/t = 0.07, where the system is always in the symmetric regime. 



Tc. In particular, the fact that local antiferromagnetic order is 
found at higher temperatures than superconducting order for 
0.6 < \/j,\/t < 0.79 does not imply that the system shows anti- 
ferromagnetic long-range order for these values of fi. For low 
enough temperatures superconducting order actually prevails 
for ji/t > 0.66 The flow equations for the regimes with spon- 
taneous symmetry breaking are discussed in more detail in the 
following section. 
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FIG. 5: Same as Fig. H now for U/t = 3, f/t = -0.1, ji/t = -0.83 
and T/t = 0.011. 



The flow of the bosonic mass terms and Yukawa couplings 
in the symmetric regime is shown in the first and second pan- 
els of Fig. m Here we have chosen a chemical potential where 
antiferromagnetism is the dominant instability. Since the ^- 
wave superconducting mass term falls slightly below the an- 
tiferromagnetic mass term and the Yukawa coupling in the d- 
wave channel hd rises above the Yukawa coupling in the anti- 
ferromagnetic channel, one has to look at the ratios /fh^ in 
order to see that the coupling in the antiferromagnetic chan- 
nel is actually the dominant one. This is shown in the third 
panel of Fig. (H where one can see that for the given choice 
of parameters the antiferromagnetic coupling is more strongly 
enhanced than the couplings in the s- and J-wave supercon- 
ducting channels. The coupling in the charge density channel 
grows least of all four. We observe the very small value of the 
eff'ective coupling for J- wave superconductivity at short dis- 
tance scales (large k). This reflects the fact that this coupling 
is only generated by the antiferromagnetic fluctuations. 
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In Fig. [5] the flow of the bo sonic mass terms, Yukawa 
coupHngs and efl'ective fermionic four-point coupHngs is dis- 
played for a combination of parameters where the coupHng in 
the J- wave superconducting channel is the dominant one. Al- 
though this coupling is smallest on high scales of the flow by 
several orders of magnitude, it is strongly enhanced during the 
flow due to antiferromagnetic fluctuations, as discussed within 
the context of the present framework in [41]. At temperatures 
slightly lower than in Fig. [5] the mass term reaches zero 
and the J-wave coupling diverges at a nonzero renormaliza- 
tion scale k = ^ssb • 




)0<>:<XX>::X 
XX' )o( >c< 



B. Quartic bosonic couplings 

The flow of the quartic bosonic couplings Aa, Ad, and Aad 
is crucial for the long-range physics of the system in the 
symmetry-broken regimes, which is dominated by bosonic 
fluctuations. In order to obtain at ^ = kssB the appro- 
priate starting values for the flow of these couplings in the 
symmetry-broken regimes, however, one has to consider their 
flow already in the symmetric regime. If commensurate anti- 
ferromagnetic fluctuations dominate, the flow equation for the 
antiferromagnetic quartic coupling Aa is given by 



p 



(p^(P)p^(p + n)) 



(30) 
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where AF^"^^ denotes the one-loop contribution to the bosonic 
four-point function, obtained as the fourth functional deriva- 
tive of the flowing action with respect to the field a, and the 
dot ■ indicates the insertion of kdk under the measure of the 
loop integral implicit in AF^'^^ Where incommensurate fluc- 
tuations dominate over commensurate ones the flow equation 
([30b for Aa has to be modified, yielding 

kdk'Aa = \ (Af -2,, 2,, -Q,) 

+Afl4)(2,,-2,,2^,-2^)). (31) 

For the quartic coupling ~Ad of the J-boson one has the flow 
equation 



kd{Ad^\Qi%^)Y/^^ 



p (4(P)P^(-P)) 



(32) 



'4 



3 

-\- -- 



The flow equation for the quartic coupling Aad describing 



FIG. 6: Contributions to the flow of the quartic bosonic couplings 
Aa, Ad and Aad in SYM. The first line shows the contributions from 
fermionic, the second and third lines those from bosonic loops. Wig- 
gly lines denote antiferromagnetic, dashed lines superconducting 
propagators, so the diagrams with four external wiggly lines con- 
tribute to the flow of Aa, those with four external dashed lines con- 
tribute to the flow of Ad and those with two external wiggly and two 
external dashed lines to the flow of Aad- 



the mutual interaction between the a- and J-boson is given by 



-^fdivf 



{p'^(P)) p^(-p)4(p + n) 



fd(v)fd(V + tt) 



pUp)pU-p)pUp + i])pU-p + n), 



AaA 



ad 



+ 2- 



ad 



{p'a(P) + ml) {p'^(P) + fhl) 



+2 



^Id 



{pi(P)^ml){p'^/p)^ml) 



(33) 



Graphical representations of the diagrams from which the 
contributions to the flow of Aa, Ad and Aad are obtained as 
scale derivatives are given in Fig. [6l 

The upper panel of Fig. [7] shows the flow of the quartic 
bosonic couplings Aa, Ad and Aad for the same set of param- 
eters as used in Fig. HI Although antiferromagnetism is the 
dominant instability for this choice of parameters, the quar- 
tic coupling Aa (green, dashed curve) is only comparatively 
weakly enhanced during the flow. For smaller values of -f 
and not so close to half filling it may even turn negative dur- 
ing the flow so that the efl'ective potential, according to the 
truncation (|23]) , is no longer bounded from below so that the 
truncation is no longer adequate and has to be replaced by a 
more extended one. A negative value of Aa may either indicate 
a tendency toward a first order antiferromagnetic phase transi- 
tion, but it may also result from a more general inadequacy of 
the parametrization of the eff'ective potential as a polynomial 
in the fields in the given range of parameters. To avoid these 
difl&culties, which do not arise at larger values of -yu (see the 
green, dashed curve in the lower panel of Fig. [7]), we focus in 
this work on values of the parameters f and ji for which Aa is 
non-negative on all scales. 
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FIG. 7: Upper panel: Flow of the (unrenormalized) quartic bosonic 
couplings ~Xa (green, dashed), Xd (red, solid) and Xad (blue, dotted) in 
SYM for Ult = 3, f It = -0.1, jilt = -0.6 and T/t = 0.07. Lower 
panel: The same fox jilt - -0.83 and T It - 0.01 1, with multiplied 
by 0.1. 



While the coupling Aa stays rather small during the flow 
and mostly has only a mild influence on the flow of the anti- 
ferromagnetic mass term in SYM, the quartic coupling Xd can 
grow very large. Already for the parameters used in the upper 
panel of Fig. [71 where the J- wave channel is far from crit- 
ical, the coupling Xd (red, solid curve) is substantially more 
enhanced than the quartic coupling Aa. The increase of Ad is 
even stronger in the range of parameters where J- wave super- 
conductivity is the dominant instability. This is shown in the 
lower panel of Fig. [71 where Ad is displayed after division by 
ten. The eminent growth of Ad during the renormalization flow 
is chiefly responsible for the fact that the transition to J- wave 
superconductivity occurs only at rather large values of -jd as 
compared to the results in Ref. '47 where no quartic bosonic 
couplings were taken into account. 



The quartic coupling Aad^ which describes the direct inter- 
action between the a- and d- boson can change its sign from 
positive to negative, or inversely, during the renormalization 
flow, see Fig. [7] (short-dashed curves). If it is positive, it en- 
hances the mass terms and m^, otherwise it decreases them 
like the fermionic contributions to their flow. 



C. Anomalous dimensions and wave function renormalization 

For the long distance behavior of the system, the anomalous 
dimensions ria and rid are of importance. They are defined as 

ria = -kdk In Aa m = -kdk In , (34) 

so they can be determined from the flow equations for A« and 
Ad. A description of how we access these quantities in the 
present approach can be found in Appendix B of Ref. 41. 

The flow equation for the fermionic wave function renor- 
malization factor Zf = Zf(oj = ±7iT) is obtained from the 
flow of the fermionic propagator at the lowest two Matsubara 
modes ±7iT . We use the formula 

kdkZF = ^ (Af^'^(7rr, q^) - Atf(-nT, q^)) . (35) 

Here the subscript F and the superscript in AF^^ indicate 
that the derivative has to be taken two times with respect to 
the fermionic fields. Again, the dot ' indicates the insertion of 
kdk under the measure of the loop integral implicit in AF^\ In 
our ansatz some choice has to be made for the Fermi momen- 
tum qf appearing on the right hand side of Eq. ([35]) . As we 
have checked, the increase of Zp during the flow is in general 
stronger for close to the points (0, ±7r) and (±7r, 0) than for 
qF close to the Brillouin zone diagonal |57], but the precise 
choice does not matter for the semi-quantitative features of 
the phase diagram. For the results displayed in the figures we 
have set qp = (0,7r). 

The flow of the Z- and A-factors used in the parametrization 
of the a- and J-boson propagators is displayed in the upper 
panels of Fig. [H The lower panel shows the fermionic wave 
function renormalization factor ZpijiT), which start its flow 
from 1 and grows by some fraction for which the increase by 
20% in Fig. [His representative. 

V. FUNCTIONAL RENORMALIZATION FOR THE 
SYMMETRY BROKEN REGIMES 

In the studied parameter region we observe spontaneous 
symmetry breaking only in the antiferromagnetic and J- wave 
superconducting channels. For the range of k where one of 
these channels shows local order we drop the charge density 
and ^--wave superconducting bosons from our truncation and 
neglect the fermionic self-energy corrections. Furthermore, 
we restrict our attention to temperatures T > Tmin = 4 • I0~^t 
and thus do not assess the ground state properties of the model 
at difl'erent values of ju and f. Furthermore, we neglect the 
scale- and momentum-dependences of the Yukawa couplings 
ha and hd, keeping their values fixed at those which they have 

at kssB- ha\k = ha(0)\kssB and hd\k = hd(0)\kssB ^r k < kssB- 
This neglect is made mainly due to computational reasons, 
but it should not have an important impact on the flow in the 
S SB -regimes at low temperatures, which is dominated by the 
long-range bosonic fluctuations, i. e. the bosonic masses, or- 
der parameters and quartic couplings. 

As a final simplification, we neglect the incommensurabil- 
ity q in the S SB -regimes, which would otherwise have to be 
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FIG. 8: Upper panel: Flow of the bosonic wave function renormal- 
ization factors (green, dashed) and Zj (red, solid). Middle panel: 
Flow of the gradient coefficients A« (green, dashed) and Aj (red, 
solid). Lower panel: Flow of the fermionic wave function renor- 
malization factor ZpinT). All curves are for the symmetric regime at 
Ult = 3, f = -0.1, lilt = -0.6 and T/t = 0.07. 



included in the truncation (l24l) for the effective potential. Al- 
though including the incommensurability may have an effect 
on the flow of the antiferromagnetic order parameter at in- 
termediate scales, we expect that it would not influence its 
flow at low scales, which is mainly determined by the number 
of Goldstone bosons. In addition to the continuous symme- 
try associated to the antiferromagnetic order parameter, in- 
commensurate antiferromagnetic order breaks the symmetry 
of rotations of the lattice by 7r/2. Nevertheless, the spon- 
taneous breakdown of this discrete symmetry does not lead 
to the emergence of additional Goldstone modes. We also 
do not think that the spontaneous breaking of lattice trans- 
lation invariance results in major changes of the flow. Set- 
ting the incommensurability to zero in the spontaneously bro- 
ken regimes will presumably leave the universal aspects of the 
flow of the running couplings in these regimes intact. 



A. Flowing potential for bosons 

For the S SB -regimes we derive the flow equations for the 
order parameters and quartic couplings from the flow equation 
of the local effective potential Ub{(^, ^) which is given by 



(36) 



Here Pk[o^, S] corresponds to the cutoff-dependent full inverse 
propagator [a, S] -\-Rk evaluated for constant bosonic fields. 
Within our truncation, the right hand side of the flow equa- 
tion for the effective potential ([36b can be decomposed into a 
fermionic and a bosonic contribution, 

dkUBia, 6) = (dkUsia, S)f + (dkUB(a, S)f . (37) 

The bosonic part can be written as 

(dkUsla, S]f = i y In [Pi(P)Sij + Mfj(a, S) + (P)^,-,,] 

^ P,iJ 

(38) 

where Pi(P) = Pa(P) and Ri(P) = Ra(P) for / = 1, 2, 3, and 
Pi(P) = p^(P) and Ri(P) = Rd(P) for / = 4,5, respectively. 
The matrix Mfj(a, S), which has to be diagonalized, has en- 
tries 



Mlj(a,S) = 



Aa(a - ao) + Aad(S - So), 
Ad(3S - So) + Aadio! - o^o), 
Ad(S - So) + Aadio^ - o^o), 

^Aad V^, 
jAad V^, 

0, 



= j=U 
-J = 2,3, 
= j = 4, 

= 7 = 5, (39) 
-l&j = 4, 
= 4 & J = 1 , 
otherwise . 



The first and fourth lines and columns of the matrix Mfj(a, S) 
are associated to the radial, the others to the Goldstone modes. 
The radial modes of the two bosons are coupled to each other 
through the coupling Aad whereas the Goldstone modes re- 
main unaffected. The form Eq. ([39b for the matrix Mf.(a, S) 
is adequate only in SSBad where the minimum of the effective 
potential U(a, S) occurs at nonzero values ao. So of both order 
parameters a and S. For ao = 0, the first three diagonal entries 
of Mfj(a, S) have to be replaced by fhl -h 3Aaa (for / = 1) and 
ml -h Aaa (for / = 2, 3). For = 0. the fourth and fifth diago- 
nal entries of Mfj(a, S) have to be replaced by ffi^ -h 3AdS and 



ml + AdS. 



The fermionic part (dkUsia, S))^ of the flow of the effective 
potential is given by 



(dkUB)' 



--Trp^kln^, 



(40) 



where the sum in the trace Trp is over fermionic indices 
only. Introducing the antiferromagnetic gap = 2hlao 
and the (momentum-dependent) J- wave superconducting gap 
^^(q) = 4/i^/^(q)^^o this fermionic contributions to the flow 
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of the effective potential can be derived from 
(AUef = -^TrFlnf' 



(41) 



= -rr^yincoshf^), 



where 
0. 



-(^p + ^p+^) + e y -(^p - ^p+^)2 + A2 



+A^(q) 



1/2 



(42) 



By the help of Eqs. ( l38t and ( l40b the flow equations for the 
quartic couplings Aa, Ad and Aad are obtained by appropriate 
derivatives with respect to the fields a and 6 on both sides of 
Eq. dlSl, 



c^/fc^aJ = — idkU(a,6))\ , ^ 



(43) 



These formulas are also valid if one of the symmetries remains 
unbroken in which case one has to set either ao or to zero. 
For ao = or So = the flow of the mass term obeys 



da 



(dkU(a,S)) 



or 



dkfh', = —(dkU(a,S))l__ 



(44) 



(45) 



While in the symmetric regime the fermionic contributions 
to the flow of the mass terms are always negative and drive the 
masses toward zero, in the regimes with spontaneous sym- 
metry breaking they may change sign. In this case even the 
fermionic contribution can lead to an increase of the mass 
terms of the bosonic fields with vanishing order parameters. In 
particular, if the antiferromagnetic order parameter acquires a 
nonzero value, this may change the sign of the fermionic con- 
tribution to the flow of the superconducting mass term and 
prevent it from becoming zero. This efl'ect is shown in Fig. 
[21 In that sense the presence of antiferromagnetic order in 
the system has a tendency to prevent the establishment of d- 
wave superconducting order. Similarly, in the regimes where 
both ao and are nonzero, a large value of ao has a dimin- 
ishing influence on the fermionic contribution to the flow of 
^0, which therefore grows less quickly or decreases faster for 
^ ^ than if ao were zero. This efl'ect acts against the co- 
existence of antiferromagnetic and J-wave superconducting 
order. Indeed, the phase diagram in Fig. [T] shows no region of 
coexistence of both orders, in contrast to what one might have 
expected from the flow of the masses and Yukawa couplings 
in the symmetric regime and the pseudocritical temperatures. 
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FIG. 9: Flow of the J- wave superconducting mass term (red, 
solid) and (unrenormalized) antiferromagnetic order parameter ao 
(green, dashed), the latter multiplied by ten. Here nonzero in- 
verts the sign of the fermionic contribution to the flow of so that 
it no longer decreases but rather increases and later saturates during 
the flow. Parameters chosen are Ult - 3, t' It - -0.1, \ilt - -0.65 
and Tjt = 0.04 



B. Flowing minimum 

To derive the flow equations of the order parameters ao and 
^0 we use the condition that Ub{ol^->So) should be the mini- 
mum of the efl'ective potential t/fi((2, &). From the necessary 
condition 

^^t/fi(ao, ^o) = dsUBi^cco. ^o) = , (46) 
which has to hold on all scales, one obtains the prescription 



d d 
—daUB(ao,So) = —dsUB(ao,So) = 0. 
dk dk 



(47) 



Together with Eq. (1431) the flow equations for the order pa- 
rameters follow: 



dkao 



dk^o = - 



^ad 



^Id 



^a^d ■ 



d,dkUB,k{a,6)l^^^^^^^^ 
d,dkUB,k{a.S)l^^^^^^^^, 
dsdkUB,k{a,S)l^^^^^^,^ 
dAUB,k{a.S)l^^^^^^,^^ 



(48) 



For parameter regions where AaXd - reaches zero the poly- 
nomial approximation for the flowing potential Ub((^, S) is no 
longer appropriate, and we discuss this issue below. 

For the studied temperature regime T > Tmin the lowest 
Matsubara mode ^ = dominates in the spontaneously bro- 
ken regimes (k < ksss) and the dimensionality of the problem 
is efl'ectively reduced from 2 -h 1 to 2, a mechanism which is 
known as "dimensional reduction". Within our flow equation 
approach dimensional reduction occurs automatically and in 
a smooth way due to the efl'ective form of the flow equations 
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1581] . It occurs effectively already for k somewhat above kgsB 
where the contribution of bosons with n becomes small. 
For computational simplicity, we therefore neglect contribu- 
tions from all bosonic Matsubara modes except the lowest 
ones in the spontaneously broken regimes. This assumption 
becomes exact in the limit of ^ «: nT, and we expect it to 
involve only a small quantitative inaccuracy at scales close to 
the critical scale. For very low temperature T < kssB these 
arguments are no longer valid. This is the reason why the 
phase diagram in Fig. [T]is not shown for temperatures close 
to r = 0. If ao or but not both, are nonzero, the long- 
range behavior of the system at finite temperature can be de- 
scribed by the 0(3)-symmetric |60] or 0(2)-symmetric lin- 
ear cr-model, depending on whether ao or is nonzero. The 
properties of these models are well-known and well under- 
stood. 

In a regime where the lowest Matsubara mode dominates 
over the others by far, one is dealing with an eff'ectively two- 
dimensional problem, the unrenormalized field expectation 
values ao and have to vanish in the infrared limit ^ ^ 
in accordance with the Mermin- Wagner theorem. However, 
for the 0(2)-symmetric model, which can be used as an ap- 
proximation in the dimensionally reduced regime for super- 
conducting order, the renormalized field expectation value 
Kd = t^AdSo/T may remain nonzero even if drops to zero as 
the gradient coefficient may diverge in this case |59]. This 
behavior is characteristic of a Kosterlitz-Thouless phase tran- 
sition [61], for functional renormalization group treatments 
see 1 38, 62, 63]. Although the polynomial expansion of the 
effective potential in Eq. (l24l) is not sufficient to account for 
the finiteness of Kd down to ^ = 0, it is sufficiently accurate 
to describe its being nonzero down to scales k ^ k^h much 
smaller than any realistic inverse probe size /"^ see 11381] for a 
more detailed discussion. 

In order to make contact with the familiar results from the 
(9(A^)-symmetric models, it is convenient to introduce the di- 
mensionless (renormalized) quantities or, /c^, 5, Kd, Aa, ^d and 
Xad, which are given by 



a = 



6 



T 

t^Ad 
T 
T 



a, 



t'Aa 

t^Ad^ 

Kd = -J^^O, 



fk^Al 
T 



Ad = 



fik^AgAd 



(49) 



from the 0{2)- and 0(3)-symmetric linear cr-models, namely 
, . (4-^J / 3 \ 

kdkKa = —- — ^ ^l\-^aKa. (50) 

167r \(\^lX^Kaf I 



kdkAa = 



^71 \{\+2AaKa? 
-2(1 - ria)^a 



+ 2 



(51) 



for the a-boson, and 



kdkKd 
kdkAd 



= A 



C^ -rid) 

167r 1(1+2^^^^)^ 
2(4-^j), 



-^l]-r]dKd, (52) 



87r \(l^2AdKd)' 
2(1 - r]d)Ad 



-h 1 



(53) 



for the J-boson. Since in the regime with two nonzero or- 
der parameters the absolute value of Aad is normally driven to 
zero much faster than the two other quartic couplings Aa and 
Ad, the flow of Ka, Kd, Aa, Ad is generally well described by 
Eqs. ([5Qb-([53l). If, however, \Aad\ is larger than the geometric 
mean of Aa and Ad, that is if \Aad\ > ^Aa • Ad, the efl'ective 
potential U (a, S) no longer has a minimum at (ao, ^o) and this 
signals the breakdown of our truncation which relies on an 
expansion of U(a,S) around (ao,So), assumed to be the lo- 
cation of a minimum. Fortunately, however, our numerical 
results for the truncation Eq. ([24l) yield a violation of the con- 
dition \Aad\ < ^Aa ' Ad only in regions where antiferromag- 
netism strongly dominates over J- wave superconductivity. In 
this regime, the effect of J- wave superconducting ffuctuations 
on the emergence of antiferromagnetic order is negligible and 
the truncation Eq. (l24l) is not natural. Consequently, if in this 
regime \Aad\ rises above ^JAa • Ad, we set Aad to zero on all 
scales whereby the expansion for the effective potential be- 
comes again well-defined. 

The main difference between the flow equations for Ka and 
Aa on the one hand and Kd and Ad on the other concerns the 
"-f2" in Eqs. (l50l) and (|5T1) as opposed to the "-Fl" in Eqs. (l52b 
and (1531) . This corresponds to the different numbers 2 and 1 of 
Goldstone bosons in the symmetry broken phases of the 0(3)- 
and (9(2)-symmetric linear cr-models, respectively. Since in 
the presence of a non-negligible order parameter the Gold- 
stone modes have a much stronger influence than the radial 
modes in driving the order parameter to zero, their number 
is crucial for how long (in terms of the renormalization group 
flow) the system remains in the symmetry broken regime. The 
beta-functions for Ka and Kd are qualitatively different, since 
for Kd the contribution "-hi" is canceled by the anomalous di- 
mension, as we will see next. 



C. Anomalous dimensions 



In terms of these quantities, when the fermions are fully 
gapped and dimensional reduction is efficient so that k <^ 2nT 
and k ^ 71, the flow equations for the order parameters and 
quartic couplings at vanishing Aad = reduce to those familiar 



In order to obtain the anomalous dimensions, one has to 
determine the flow equations for Aa and Ad in the presence of 
nonzero Ka and/or Kd. To this end, we take a second derivative 
of the loop contributions to Pa(Q) and Pd(Q) with respect to 
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spatial momentum and a derivative with respect to the scale k: 

q2 



^,A, = ^,|limi^AP,(0,/,0)j, 



(54) 
(55) 



In the regimes exhibiting spontaneous symmetry break- 
ing, the fermionic contributions to t/^ and rid quickly become 
negligible as soon as the scale drops below the temperature, 
and it suffices to consider the bosonic contributions. In case 
the two bosons can independently be described by the two- 
dimensional 0(3)- and 0(2)-symmetric models these contri- 
butions are, assuming dimensional reduction. 



1 

na,d = - 



2 ' 



(56) 



In the presence of nonzero Aad, this formula has to be gener- 
alized, yielding 



f<d(^ - ^f<a(Aa - KdAl, + iKdAaAd)) 



ad 



(l -h 2KdAd + 2Ka(Aa - ^KdA^^ + iKdAaAd)) 
KaAl(l + iKdAdf 



(l + iKdAd + 2Ka{Aa - + 2KdAaAd)) ^ 

^ Ka(l - 4Kd(Ad - KaA^^ + 2KaAaAd)) 



,(57) 



"^ad 



(l -h 2KdAd + 2Ka{Aa - + 2KdAaAd)) 



KdA]{l+2KaAaf 



(l + 2KdAd + 2Ka{Aa - + 2KdAaAd)) ^ 

which reduces to Eq. ([56b for Aad = 0. 



(58) 



f<a 



lilt = -0.4 
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VI. TEMPERATURE DEPENDENCE OF THE ORDER 
PARAMETERS 

We now turn to the discussion of our numerical results for 
the order parameters and the gaps as functions of temperature 
(Fig. fTOl) , the flow of the running couplings in a selected num- 
ber of cases (Figs. [TT]and[T2]), and the more general features 
of the phase diagram (Fig. [TJ in the symmetry broken regimes 
SSBa, SSBd and SSBad. 



A. Order parameters 

Fig. [TOl displays the renormalized antiferromagnetic and 
J- wave superconducting order parameters Ka and Kd as func- 
tions of temperature at difl'erent values of the chemical poten- 
tial II for U = 3t and f = -OAt. Both Ka and Kd are eval- 
uated at kph = lO'^t ^ 1 cm"\ corresponding to a realistic 
inverse probe size. The upper panel shows the temperature 



dependence of ka = Ka 



at the van Hove filling // = 4f . 



FIG. 10: Renormalized order parameters ka and kd at the "macro- 
scopic" scale k = kph, corresponding to an inverse probe size of 
^ 1 cm, as a function of temperature T for U/t = 3 and different 
values of fi. The upper panel shows the temperature dependence ka 
for II It - -0.4, the middle panel shows kd for njt - -0.81. The 
lower panel shows kd and ka for iijt - -0.72, where they are nonzero 
in different temperature ranges. 



The shape of this curve for ka is similar to that of the curve 
presented Fig. 1 in Ref. 36 for = n = 0. The temper- 
atures where ka is nonzero, however, are lower according to 
the results presented here since more fluctuations have been 
included which have a tendency to destroy antiferromagnetic 
order. The middle panel of Fig. [TOl shows the temperature 
dependence of kd = KdV. dX ult = -0.81, where only d- 

IK—Kph 

wave superconducting order and no antiferromagnetic order 
occurs. The steep fall to zero of kd Sit T = Tc can be seen 
as a remnant of the jump in the superfluid density found for 
a Kosterlitz-Thouless phase transition at Tc in an improved 
truncation [38^ .631 
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In the absence of electromagnetic interactions (as for the 
pure Hubbard model) the superconducting phase is actually a 
superfluid phase with superfluid density given by Us = Kdla^, 
where a denotes the lattice spacing and Us the number of par- 
ticles per area. For nonzero electromagnetic coupling e the 
"photon mass" responsible for superconductivity is given by 
My = 2e^fkd. (Here we observe that the J-boson carries 
charge two. More precisely, e is the effective renormalized 
electromagnetic coupling at the scale kph.) Other observ- 
able quantities are the effective gaps for the electrons. In- 
deed, the antiferromagnetic and J- wave superconducting gaps 
and are related to ka and kd by A^ = ^|2h^ and 

Aj(q) = 'yjh^fd((i)^kd, where the renormalized Yukawa cou- 
plings ha and hd are given by /z^^ = Ah^^l,d' ^^S- IIHl there- 
fore predicts measurable quantities. 

The lower panel of Fig. [TOlshows a situation where, Sit/u/t = 
-0.72, we find nonzero Kd at low temperatures and nonzero 
Ka at higher temperatures. In between, we observe a small 
temperature region around T = 0.0 125^ where neither of the 
two order parameters remains nonzero down to ^ = kph. In 
this region only local but no long-range order is present. An 
interesting feature of this graph is the steepness of the rise 
of ka 3.1 T = Tc, which contrasts with the behavior at van 
Hove filling (upper panel of Fig. [TOb where the rise below 
Tc is relatively smooth. The main reason for this feature is 
the strong growth with temperature of the final value of Aa in 
SYM close to r = Tc. This value has an important influence 
on the initial growth of Ka in the spontaneously broken regime 
and therefore at its value at ^ = kph. This effect is mainly 
responsible for the smallness of the temperature interval in 
which ka drops to zero as the temperature approaches Tc from 
below. 

The fact that both order parameters become zero during 
the flow at values of k > kph in a temperature region around 
T = 0.01 25 ^ results from the mutual negative influence of 
the two types of order on each other. This influence is fur- 
ther illustrated by the upper and lower panels of Fig. [TTl 
where the upper panel shows the flow of Ka together with that 
of Kd down to ^ = kph. For the temperature T = 0.0131^ 
used in this graph Ka is nonzero during a much longer pe- 
riod of the renormalization flow than /cj, which becomes zero 
at - ln(^/0 ^ 4.6 so that the superconducting mass term ffi^ 
(magenta curve) becomes nonzero again. The lower panel of 
Fig. [m in contrast, shows the flow of the same couplings, but 
in this case the mutual influence of the order parameters has 
been neglected. This means that each order parameter is set 
to zero in all contributions to the other boson and the inter- 
bosonic quartic coupling Aad is set to zero. As described in 
the previous section, this is equivalent to deriving the bosonic 
contributions to the flow equations from the 0(3)- and 0(2)- 
symmetric linear cr-models at finite temperature. According 
to this simplified treatment, neglecting the mutual influence of 
the two types of order, both order parameters remain nonzero 
down to ^ = kph. Such a result would suggest a region of 
coexistence of "global" antiferromagnetic and J- wave super- 
conducting order. 

In our example this coexistence is destroyed by the mu- 
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FIG. 11: Upper panel: Flow of the renormalized antiferromagnetic 
order parameter Ka (green, dashed), renormalized J- wave supercon- 
ducting order parameter /cj (red, solid) at U/t = 3, f /t = -0.1, 
jd/t = -0.72 and T/t = 0.0131. The dashed-dotted magenta curve 
shows the (unrenormalized) J-wave superconducting mass term rfi^ 
when it becomes nonzero again. Lower panel: Same as upper panel, 
but neglecting the mutual influence of the order parameters, i. e. each 
order parameter is set to zero in all contributions to the other boson 
as well as the inter-boson coupling Aad. 



tual influence of the antiferromagnetic and superconducting 
bosons. It is precisely this type of influence which has been 
taken into account in the upper panel of Fig. [TTl We there- 
fore conclude that the two types of order have a tendency to 
destroy each other. On the basis of a renormalized mean field 
treatment Ref. (see in particular Figs. 10 and 11) reports 
on an analogous tendency of antiferromagnetism and super- 
conductivity to mutually suppress each other. 

For the curves shown in Fig. [121 the temperature has been 
reduced in comparison to Fig. [TTl so that both order param- 
eters (upper panel) are nonzero for an important interval of 
the flow. Now J- wave superconducting order persists down 
to much lower scales k of the renormalization flow. Although 
at an intermediate stage of the flow Ka is considerably larger 
than Kd, it vanishes earlier during the flow due to the larger 
number of Goldstone modes for antiferromagnetism. In the 
lower panel of Fig. [121 the flow of the quartic couplings Aa, Ad 
and Aad is displayed, where Aad approaches zero much more 
quickly than Aa and Ad so that the two bosons are more or 
less independent and the flow is dominated by their Goldstone 
modes at low scales. 

Taking things together, the example shown in the lower 
panel of Fig. [TOl demonstrates that the phase transitions can- 
not always be understood by the universal behavior of linear 
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FIG. 12: Upper panel: Flow of the renormalized antiferromagnetic 
order parameter Ka (green, dashed) and renormalized J- wave super- 
conducting order parameter /cj (red, solid) dXUIt - 3, f It - -0.1, 
jilt - -0.72 and T It - 0.0118. Lower panel: Flow of the quartic 
couplings Aa (green, dashed). Ad (red, solid) and Aad (blue, dotted), 
the latter two multiplied by a factor of 0.1, for the same choice of 
parameters. 



or non-linear uncoupled cr-models. For example, in the (9(3)- 
(T-model it is not possible to find a restoration of disorder at 
temperatures below the ones for which long-range order is re- 
alized. The competition of diff'erent bosons is crucial for a 
quantitative understanding of the phase diagram. 



B. Phase diagram 

We now turn to the discussion of the phase diagram ob- 
tained for U = 3t and f = -0.lt, as shown in Fig. [T] For 
values of -f which are substantially smaller than 0.1^ the 
quartic coupling Aa eventually becomes negative during the 
flow, which may indicate a tendency toward a first order tran- 
sition which is not captured in the present truncation of the 
efl'ective potential (l24l) . For values of -f which are consider- 
ably larger, in contrast, the system exhibits a tendency toward 
ferromagnetism (y\\\ . In order to account for this instability, 
the truncation for the efl'ective action and the parametrization 
of the bosonic propagators and Yukawa couplings specified 
in Appendix B of ll4lll would have to be adjusted accord- 
ingly. Upon small variations of f the qualitative picture of 
the phase diagram remains essentially unchanged. If -f is re- 
duced, all phase boundaries are shifted toward smaller values 
of -yu, if -f is increased, they move into the other direction. 



For smaller values of the Hubbard interaction U, critical tem- 
peratures are lower and the phase boundaries are shifted in the 
direction of the van Hove filling chemical potential ji = 4t\ 
The results we have obtained for calculations at values of U 
and f other than U = 3t and f = -0.lt do not alter the picture 
described in what follows. 



At the van Hove filling we find a sizable diff'erence be- 
tween the pseudocritical temperature Tpc and the true crit- 
ical temperature Tc for antiferromagnetism, which diff'er by 
a factor of about 2, mainly due to the two antiferromagnetic 
Goldstone modes. In the J- wave superconducting regime at 
-j^lt > 0.75, in contrast, there is only a slight diff'erence be- 
tween Tpc and in accordance with earlier results on the 
0(2)-symmetric model |38]. The non-negligible diff'erence 
between Tpc and Tc for J- wave superconductivity in the region 
between -j^lt = 0.66 and -j^lt = 0.75 is not due to Goldstone 
fluctuations but arises from the influence of antiferromagnetic 
order on the flow of /cj. 



One of the most intriguing questions about the phase di- 
agram of the two-dimensional Hubbard model is whether it 
exhibits a region in parameter space where antiferromagnetic 
and J-wave superconducting order coexist. In principle the 
setup employed in the present work makes it possible to as- 
sess this question, but the results obtained by means of the 
truncation used here do not permit a definite answer. While 
they clearly suggest that the two types of order "do not like 
each other", they can hardly be taken to rule out the existence 
of a region in parameter space where antiferromagnetism and 
J- wave superconductivity coexist. Where the two phases bor- 
der each other at low temperatures around -fi/t = 0.66 in 
Fig. [T] there is always at least one type of order which persists 
down to k = kph = I cm"^ ^ 10"^t, and the value of k where 
the second order parameter drops to zero is often very close 
to kph. Coexistence might occur, even within the truncation 
used here at temperatures below Tmin = 4 • 10"^^, the lowest 
temperature for which we have done calculations. In all cases 
where both order parameters remain finite for a considerable 
part of the flow, the values of the running couplings at ^ = kph 
are highly sensitive to their values at the onset of the spon- 
taneously broken regime. Therefore, we expect that further 
extensions of the truncation, which may influence the flow on 
intermediate scales, can have an important efl'ect on the shape 
of the phase boundaries where the antiferromagnetic and d- 
wave superconducting phases are close to each other. Self- 
energy corrections, higher order bosonic couplings and the 
efl'ect of the antiferromagnetic incommensurability, which we 
have neglected here in the SSB regimes, may be responsible 
for whether there exists a region in the phase diagram where 
antiferromagnetic and J- wave superconducting order coexist 
at ^ = kph. The results presented here, however, suggest that 
if there is a region in ji - T-space where antiferromagnetism 
and J- wave superconductivity coexist on a macroscopic level, 
this region is probably not very extended. 
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VII. CONCLUSIONS 

Our main conclusion is that functional renormalization 
combined with partial bosonization can give direct and phys- 
ically transparent access to the low temperature behavior of 
the Hubbard model. We have computed the temperature de- 
pendence of the superfluid density or the antiferromagnetic 
order parameter for temperatures below the critical temper- 
ature from which the associated gaps for the electrons can 
easily be determined. Taking into account the fluctuations of 
composite bosons we can incorporate the important contribu- 
tions from collective spin waves or J-wave superconducting 
bosons. This allows us to compute the critical temperature as 
a function of the chemical potential and therefore to establish 
the phase diagram. 

The physics of collective boson fluctuations is quantita- 
tively important. This is demonstrated in Fig. [T] by the 
substantial difl'erence between the pseudocritical temperature 
(which is often associated with the critical temperature in 
other approaches) and the true critical temperature. Also the 
understanding of the region with coexisting local antiferro- 
magnetic and superconducting order (in Fig. [T]for -ji/t be- 
tween approximately 0.6 and 0.8) would be difl'erent without 
a proper understanding of the bosonic fluctuations. 

It is obvious that more extended truncations can substan- 



tially improve the quantitative accuracy. One may include 
into the flowing momentum-dependent four-fermion vertex 
(instead of a constant U) the changes which cannot be ab- 
sorbed by bosonization into the present channels, or one could 
increase the number of bosonic channels retained. Further- 
more, a more accurate parametrization of the momentum and 
frequency dependence of the bosonic and fermionic propaga- 
tors could be helpful. One may go beyond the quartic polyno- 
mial approximation for the efl'ective bosonic potential, both in 
order to account for possible first order phase transitions and 
to give a more accurate description of the universal critical 
behavior for the second order phase transitions. 

Nevertheless, we believe that several of our findings are 
rather robust. This concerns the existence of an extended re- 
gion in the phase diagram for which local order but no long- 
range order exists. In particular, the transition between anti- 
ferromagnetic and superconducting order is rather complex, 
with coexisting local antiferromagnetic and superconducting 
order, but also a strong tendency of exclusion of coexisting 
long-range order due to the bosonic fluctuations. 
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